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Abstract 

We consider lepton pair production via two-hadron annihilation through 
various isoscalar vector mesons within hot, baryon-free matter. This is tanta- 
mount to constructing effective form factors which we model using a vector- 
meson-dominance approach and compare with experiment. In particular, we 
consider the reactions tt/o — > e'^e~ and KK*[S92) + c.c— s- e^e~. We find that 
u;(1390) and (/)(1680) are visible in the mass spectrum for the thermal pro- 
duction rate above the tt+tt" — > e'^e~ tail and even above the vrai — > e"'"e~ 
results — both of which were considered important in their respective mass 
regions. 

PACS numbers: 25.75.-Fr, 12.38.Mh, 13.75.Lb 
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I. INTRODUCTION 



Production of dileptons in high-energy heavy-ion coUisions continues to be studied ac- 
tively since it is hoped that such signals might be useful in discriminating prehadronic or 
quantum chromodynamic (QCD) plasma formation from ordinary hadronic matter. The pos- 
sible success depends on having a reasonably complete picture of both plasma and hadronic 
phases. Crossover from one to the other is expected for a critical temperature in the neigh- 
borhood of 200 MeV and near this boundary the hadronic matter will be populated by many 
species of strange and nonstrange mesons. However, to get a rough idea about this phase of 
matter one typically used a simple pion-gas approximation. In terms of dilepton spectra, it 
has recently been shown that thermal dielectron production rates from the pion gas could 
be an order of magnitude lower than those from a hadron gas consisting of pions and various 
resonances IQ. Going beyond simple pion gas considerations is important. Allowing heavier 
species to partially comprise the system complicates matters since each species might inter- 
act among themselves as well as with others. The list of possible processes then becomes 
quite lengthy. Further complications or uncertainties arise since, with the exception of the 
pion 0, hadronic form factors are not particularly well known. Reactions which involve 
an intermediate isoscalar vector meson could contribute to dilepton production and for this 
reason the mass distribution of produced dileptons is expected to exhibit structure owing to 
these hadronic resonances. 

Observed spectra from heavy-ion collisions will of course be restricted to a time- integrated 
yield. Temperature (or time) dependent rates are the essential quantity that might be input 
for models of the time evolution in order to make predictions for overall yields. Rate cal- 
culations have been done with various levels of sophistication. In a quest for more accurate 
knowledge of dielectron production rates we consider some hadronic processes which produce 
e+e^ pairs and then compare with spectra calculated previously The particular pro- 

cesses we consider are pseudoscalar and vector reactions of the type P + V ^ V^^^ e~^e~ 
where P and V are either pions and p-mesons or kaons and i^*(892)s. The intermediate 
isoscalar vector hadronic resonances are 0(1020), ci;(1390) and 0^(1600) for the up processes 
and 0(1680) for the strange-particle process. Special attention is payed to experimental data 
on the time-reversed processes in order to constrain or even fix the coupling constants. 

Our paper is organized in the following way. In Sect, y we present the basic formalism for 
the interactions among the hadronic degrees of freedom and the elementary scattering (an- 
nihilation) processes. A vector-meson-dominance (VDM) approach will be presented which 
allows calculation of the basic Feynman diagram for a+b — > V^^^ ^ 7* — > e~^e~ . Within this 
approach the amplitude for each diagram contains a ratio of the strong-to-electromagnetic 
coupling constants. Experimental data on cross sections for the time reversed processes 
e~^e~ ^ a + b are shown as compared with model calculations. These are used to fix cou- 
pling constants and phases between diagrams. One is really constructing an effective form 
factor by this prescription. Section |T| contains the essential details and results for thermal 
production rates which include effects from these hadronic resonances. The interesting result 
is that the uj{1390) and 0(1680) peaks are visible above the rate from vr+vr" annihilation 
and even vrai reactions which have recently been shown to be dominant for restricted masses 
0]. Finally, in Sect. |V| we conclude with some final remarks which essentially stress the 
importance of appealing to experimental data at some level for the elementary reactions in 
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future studies. Putting them into calculations on systems of hadronic matter only after such 
comparison minimizes uncertainty. 



II. FORMALISM 

After the initial stage of a collision of the type we consider, which may well include a 
prehadronic or quark-gluon-plasma (QGP) phase, hadronization gives rise to a hot ensemble 
of the lighter mesons. Within the context of a thermodynamic model this system's composi- 
tion has recently been studied . The most abundant hadronic species in this scenario are 
pions, kaons, p-mesons and if* (892) strange vector mesons owing to the large spin-isospin 
degeneracy of the last two. With temperatures ranging from phase boundary or crossover 
temperature downward toward the pion mass or lower, where the system ceases interacting, 
the number densities of these species range from several hundredths to a few tenths per fm^ 
1^. Until freezeout, they interact very strongly. Though binary reactions dominate at these 
densities, studies of the effect of ternary reactions have been pursued 0. In this work our 
aim was to study the relative strength of different lepton pair producing sources and the 
possible manifestation of some higher-lying resonances. For this, we have chosen a simple 
thermal environment on purpose. As is well known from the analysis of heavy ion collisions 
in the intermediate energy regime, nonequilibrium effects might manifest themselves in the 
final state. A related issue here is the importance of Drell-Yan and DD lepton pair signals 
around the invariant mass region adjacent to the J/V' [Q. Those are obviously important 
and should be adressed in a realistic comparison with heavy ion experimental data. 

Species can be classified according to their quantum numbers in order to identify the 
manner with which we model their interactions. Pseudoscalar particles (P) interact with a 
vector particle (V) in a manner described by the Lagrangian 

LvPP = gvPpV^Pd^P (2.1) 

and vectors interact with a pseudoscalar through 

Lvvp = 9VVP W d^Vd'^V P. (2.2) 

Let two hadrons, call them a and b, scatter in a timelike fashion, i.e. let them annihilate 
to form a resonant state V — restricted here to neutral vector mesons. It is not entirely 
clear how it converts itself from a hadronic to an electromagnetic field, but we model the 
phenomenon using the principle of vector-meson dominance. In practice, it means that the 
neutral component of the hadronic vector field couples directly to the photon via M 



em 



V 



Vl^^^A^" (2.3) 



^ 9v 

where V^^'' is the neutral component of the resonant isovector field and A is the electro- 



magnetic field. Using Eqs. ( p.l|) or (2^) and (2^) we construct the general annihilation 
diagram through a vector resonance ^ to a virtual or massive photon, which subsequently 
decays into a lepton pair as shown in Fig |l]. This differs from an analogous process of purely 
electromagnetic exchange because of the presence of the hadronic resonance. Since the effect 
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of this resonance appears as a function of invariant mass alone, the hadronic dependence of 
V at the vertex factorizes from the rest of the diagram in Fig. |l| as a form factor times a 
pointhke electromagnetic process. Vector-meson dominance has been very successful in the 
areas of hadronic form factors, photoproduction and absorption cross sections, and in the 
vector meson exchange contributions to tcN and A^A^ scattering. 

An invariant amplitude for scattering or decay in which the photon couples to the vector 
meson contains the ratio gvab/gv which forces the appearance of the same ratio (squared) 
in the cross section or decay rate. We will discuss two methods for fixing these quantities. 
For the case of the 0(1020) we can fix these coupling constants independently via strong 
and electromagnetic decays; whereas for the higher excitations presently considered we are 
forced to use another procedure. First consider the decay — > pn. We calculate its rate to 
be 
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(2.4) 



where p is the center-of-mass momentum of the decay products. This determines ^05'|p^/47r 
= 0.29 when the partial decay rate is taken to be 12.9% of 4.43 MeV. The same expression 
and numerical conclusions were reached in a recent study of 0-meson properties at zero and 
finite temperature 0. Next, the electromagnetic decay e~^e~ is modeled and computed 
to be 



(2.5) 



Using the branching fraction F^^^e+e- /r^^" = 3.09 x 10~^, we arrive at the dimensionless 
vector- dominance coupling constant for the (p of = 12.9. 

Since data are not sufficient to allow us to follow this prescription for the other resonances, 
we instead appeal to the time-reversed hadron production processes 

e+e-^a + b + ... (2.6) 

We compute the cross sections for e'^e^ pn and for e^e^ —>■ KK* {892) which is approxi- 
mated by e+e^ KgK^n^. Henceforth i^'*(892) will be written just as K*. The diagrams 
are shown in Figs. ^ and Qb. The cross section for the process in Fig. ^ is found to be 

, , 47ra2|p|3 |F^„(M)P , , 

ae+e-^.p(M) = l^lp (2-7) 



where M = \ p^ + ml. + \ p^ + m? and 



p 

with the sum running over the three vectors 0(1020), ci;(1390) and ^7(1600). ^7(783) could 
also be included in the sum but since it is so far off shell for masses considered here it can 
safely be neglected. Comparison with experiment constrains the ratios of coupling constants 
as well as phase angles in Eq. (|2.8| ). The comparison to data is shown in Fig. ^. Details 
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of the parameters in our model go as follows. For the phases we take ip^ = (without loss 
of generality one of them can be set to zero), (fi^ = dir/A and (fi^> = vr; while the ratios are 
taken to be g^j-np I Quj = 0.29 GeV~"^ and Qui'ttpIQu' = 0.05 GeV~^. Note that g(j,-Kplg<j, = 0.15 
GeV~^ was already fixed by Eqs (|2.4H 275|). Notation used here is that (j) means 0(1020), uj 
stands for ci;(1390) and uj' is ci;(1600). 

The strange particle process e^e~ — > KK* approximated by the reaction with the even- 
tual three-body final state is shown in Fig. |^. Since the decay of K* is practically entirely 
into IT + these are nearly equivalent. We take a relativistic Breit-Wigner form for the 
cross section. Computing Fg+e-^^' and Fj,/^^^,, a can be written as 



^kMM) ^^^^ \Fkk*{M)\ (2.9) 

where A(x, y, z) = — 2x{y + z) + {y — z)'^ and 



FkkAM) = [ ) — , . ^ ■ (2.10) 

Using value gcp'KK* / Qcp' = 0.19 GeV~^ we generate the curve shown in Fig. ^ compared 
with experiment. Also shown is the result from a fully microscopic calculation of the exact 
diagram of Fig. ^ whose cross section involves integration over full three-body phase space. 
Given the forms of the interactions written in Eqs. (p.l|) and ( p.2|) , the model form of the 
cross section is well defined. It is presented as the dashed curve. The same ratio of coupling 
constants is used in both approaches. 
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III. THERMAL PRODUCTION RATES 

Hadrons in this hot reaction zone are assumed to have their momenta thermally dis- 
tributed. Then the differential rate at which they scatter, i.e. the differential rate at which 
hadrons annihilate to create an e^e~ pair of invariant mass M can be written as 

dR f d^pa d^pb d^p+ d^p. 



dM^ J 2E„(27r)3 2Eb(27r)3 2E+(27r)3 2E_(27r)3' 

X (27r)^| -M \'5\pa +Pb-P+- P-)5 (m^ - (j9+ + p_f) (3.1) 

where / is the Bose-Einstein distribution and M is an overall degeneracy factor. With the 
aid of the 5 functions some of the phase space can be reduced analytically. We integrate the 
rest numerically. The absolute square of the scattering amplitude for vrp — e+e~ (initial 
spin averaged and final spin summed) is 

I r = ^ '^^ff W ^nx.r qWpP; (-/^) [plp'- + P'-Pi - 9^\P+ ■ P-)] (3.2) 

where where = p^ + p^ so that = and finally, where the form factor is given by 
Eq. (^. 

For later convenience it is useful to compute the cross section arising from this squared 
amplitude. We get 
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(T^p^e+e-(M) = ^^jj:^ ^. (3.3) 



Then the thermal rate can be written 
dR U 



dM^ 2(27r)^ 

where 



J dE^dE,f{E^)f{E,)a{s)>^''\s,ml,ml)Q{x) (3.4) 



X = ^pIpI -{ml + ml-s + 2EpE^f . (3.5) 



Notice that s and have been used interchangeably. The integration required by Eq. 
can of course be done numerically for arbitrary distributions. However, assuming Boltzmann 
distributions for the mesons it is a straigtforward exercise to perform the integration and 
arrive at the much simpler expression for the rate 

d n T 

T^TTTTJ K,{M/T)X{s, ml ml)a{s) (3.6) 



dM^ 32tt^M 

where Ki is the first modified Bessel function. Making the replacements mp — > m^* , m-j^ — *• 
mx and using the appropriate degeneracy J\f, one arrives at the thermal rate for lepton pair 
emission through KK* + c.c. e^e~ . 

Results at 150 MeV temperature presented in Fig. ^. We have compared rates which 
were integrated numerically starting from Bose-Einstein distributions to those computed 
using classical distributions written finally as Eq. ( [3.61 ). Since there is no final-state hadron 
in these reactions and therefore no Bose-enhancement effects, the classical distributions 
are more than adequate. For comparison purposes we also include the results for tt+tt" 
annihilation []1| and for vrai reactions from Ref. which were concluded to be dominant for 
masses above m^. The importance of the excitations of ci;(1390) and 0(1680) are immediate: 
they are visible above the tail of the 7r"'"7r~ annihilation and even above the vrai results. One 
last figure is shown (Fig. ^) where a temperature of 200 MeV is used. Our conclusions are 
independent of the temperature. 

IV. FINAL REMARKS 

Timelike photons are observed to have a strong affinity for neutral hadronic fields. This 
observation is captured in the current field identity [0 which is a consequence of the VMD 
hypothesis. These photons, observed as lepton pairs, exhibit resonant structure when the 
observed particles are projected onto a mass spectrum. One understands these peaks within 
vector-meson dominance to be manifestations of hadronic form factors. Our study suggests 
that in order to improve knowledge of dielectron production rates from hot hadronic matter, 
at least for dielectron masses ranging from m^ to mj/^, form factors for hadrons other than 
just the pion must be considered. We show only two examples of additional structure in the 
mass spectrum — there may well be more surprises. 

World data on electron-positron annihilation into hadrons is quite extensive at higher 
energies but is unfortunately not as complete for lower energies. However, that which does 
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exist should be exploited to its fullest potential. In our study we have constructed effective 
form factors |F^p(A'/)p and \Fkk*{M)\'^ in the mass range 1-2.5 GeV. After doing so, we 
have compared each elementary process with known data before using them in a thermal 
calculation on the extended system of interacting hadrons. Higher excitations of the lo and 
0, namely a;(1390) and 0(1680), were found to be visible in the mass distribution of the 
dilepton production rate above the hadronic production mechanisms previously considered. 
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FIGURES 



FIG. 1. General annihilation diagram of two hadrons a and 6 to a lepton pair. The process 
goes through a vector hadronic resonance V and by vector-meson dominance, goes directly to a 
virtual photon. 

FIG. 2. Annihilation of e^e^ into vrp through any of three isoscalar vectors (/)(1020), uj{1390) 
and u;(1600) in (a) and in (b) e+e" K'^K^tt^. 

FIG. 3. Cross section for e~^e~ np as computed within the model described in the text as 
compared with experimental data from Ref. |10|,11|. 



FIG. 4. Cross section for e^e~ — > KK* described by the Breit-Wigner form presented in the 
text (solid curve) , results of a microscopic calculation for the full process (dashed curve) and finally, 
experimental data for the reaction e'^e~ — > Kt;K^'iT^ from Ref. |jlz 



FIG. 5. Thermal production rate at T = 150 MeV from np (solid curve) processes, ^i^*+c.c. 
(dotted curve), vr+vr^ annihilation from Ref. |j^] (dot-dash curve) and finally, vrai process as calcu- 
lated in Ref. (dashed curve). 

FIG. 6. Same as Fig. | except using T = 200 MeV. 
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Figure 1 

K. Haglin and C. Gale, ' 'Higher Excitations of omega and phi in DileptonSpectra' ' 





Figure 2 



K. Haglin and C. Gale, ' 'Higher Excitations of omega and phi in DileptonSpectra' ' 
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K. Haglin and C. Gale, "Higher Excitations of omega and phi in DileptonSpectra" - Figure 5 
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K. Haglin and C. Gale, "Higher Excitations of omega and phi in DileptonSpectra" - Figure 6 



